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Abstract: Nonsinglet contributions to the gi(x, Q 2 ) structure function are calculated in 
the double-logarithmic approximation of perturbative QCD in the region x <C 1. Double 
logarithmic contributions of the type (a s \n 2 (l/x)) h which are not included in the GLAP 
evolution equations are shown to give a stronger rise at small-x than the extrapolation of 
the GLAP expressions. Further enhancement in the small-x region is due to non-ladder 
Feynman graphs which in the DLA of the unpolarized structure functions do not contribute. 
Compared to the conventional GLAP method (where neither the whole kinematical region 
which gives the double logs nor the non-ladder graphs are taken into account) our results 
lead to a growth at small-x which, for HERA parameters, can be larger by up to factor of 
10 or more. 
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1 Introduction 



In the framework of perturbative QCD the theoretical investigation of deep inelastic scatter- 
ing in the HERA regime puts particular emphasis on the calculation of structure functions 
in the region x<l. So far, main attention has been given to the small-x behavior of the 
gluon structure function: from the theoretical point of view, the strong rise predicted by 
perturbative QCD violates unitarity and hence requires corrections which restore unitarity. 
Experimentally, such a rise has been observed, and it is attractive to interpret it as a mani- 
festation of the BFKL Pomeron |I] . 

From the theoretical side, however, also the fermion structure functions have quite an inter- 
esting small-x behaviour. Recently it has been shown, for the case of the flavor nonsinglet 
contribution to F\, that the small-x behavior of quark structure functions is stronger than 
what one would obtain from simply extrapolating the GLAP || |J evolution equations into 
the small-x region. The reason for this lies in the fact that at small-x a new region in phase 
space opens up which gives rise to logarithms of the type (ot 3 ln 2 (l/x))™ and which is not 
taken into account within the GLAP method. As long as we are are dealing only with neu- 
tral currents and unpolarized structure functions, this observation may seem to be slightly 
academic, since the gluons in the flavor singlet will dominate at small x, and the quarks 
represent a nonleading effect. The situation, however, changes drastically if one considers, 
for example, polarized structure functions and their sum rules ||. It is well known Q that 
the polarized gluon structure function is no longer growing as (l/x) x with A close to 1, but 
rather with a A close to 0. Consequently, the polarized gluons are no longer dominating over 
the fermions, and the small-x behavior of the quarks becomes as vital for polarized deep 
inelastic structure functions as that of the gluons. 

It may be useful to briefly review the small-x behavior of the (unpolarized) gluon structure 
function in the GLAP scheme and in the BFKL Pomeron. As it is well known, GLAP sums 
logarithms of the type (a s \n(Q 2 / /j 1 2 )) k , whereas BFKL keeps terms of the form (a s ln(l/x)) fc . 
In the region of mutual overlap where the DLA is valid one is summing terms of the form 
(a s ln(Q 2 / fi 2 ) ln(l/x)) fc . In terms of the anomalous dimension 7 99 of the two-gluon operator, 
the GLAP formalism uses a fixed order (one or two-loop) in a s expression. As a function of 
n, jggin) is singular near n — 1 (n is the moment index which in the small-x region coincides 
with angular momentum j of the cross channel). In contrast, the BFKL Pomeron provides 
all singular terms of the form J2{ a s/{ n — l)) fc Cfc. Finally, in a physical gauge both GLAP 
and the BFKL Pomeron are sums of ladder diagrams, but the phase space for the transverse 
momenta of the gluons has, in the GLAP case, the property of strong ordering, whereas in 
the BFKL approximation this feature is almost completely lost Q. 

This situation changes if we consider fermions at small x 0. It has been observed al- 
most 30 years ago f| that in the Regge limit scattering amplitudes with two fermions 

3 During the rapidity evolution of the BFKL Pomeron, the mean value of the logarithm of the transverse 
momentum is still growing, but the fluctuations are much larger than in the GLAP case. 
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in the t-channel have two powers of ln(s//z 2 ) per loop, i.e. one is summing terms of the 
form (a ln 2 (s//i 2 ))™ rather than single logs. These double-logs cannot be reached within 
the GLAP scheme which, at small x, only includes the (less important) terms of the form 
(a s ln(Q 2 / /j, 2 ) ln(l/x))". This difference is seen very clearly if one analyses the relevant re- 
gion of integration inside the fermion ladders: GLAP has the same ordering of transverse 
momenta as in the gluon ladders, but the double-logs |J come from a region with no order- 
ing in kj- or the angle 9. Only the longitudinal components of the Sudakov variables are still 
ordered. In other words, in addition to the GLAP region there is another part of the phase 
space which contributes to these double-logs. Finally, it is instructive to describe this differ- 
ence also in terms of the anomalous dimension. Within the GLAP analysis the anomalous 
dimension of the quarks, as a function of the moment index n ~ j, has a pole near n=0. In 
the Regge analysis which includes the double-logs, one obtains a quite different dependence 
upon n; only for ^/oT s C n < 1 one recovers the usual result. The small-x behavior, on the 
other hand probes the region n ~ JoT s , where the two methods lead to different predictions. 




In this paper we will calculate the small-x behavior of the nonsinglet contribution to the 
polarized structure function g\\ in a forthcoming paper we will extend our analysis to the 
flavor singlet part which contains mixing between the quarks and the gluons. We find that 
the small-x behavior is stronger than the extrapolation of the GLAP formula predicts. In 
comparison with the fermions in the unpolarized structure function || the polarized quark 
structure functions have another new property. In contrast to F\ where a simple ladder 
structure gave the relevant logarithms, the polarized structure function g\ receives contribu- 
tions also from non-ladder graphs. Our analysis will show that these terms lead to a further 
enhancement at small-x. 

This paper is organized as follows. In section 2 we briefly review the definitions of the 
structure functions, and we discuss their signature properties by studying the double loga- 
rithmic contributions of the first two loop corrections. In section 3 we construct the infrared 
evolution equation for gi and find its solution. The final section 4 contains a brief discussion 
of our results. 

2 Signatures of the Structure Functions 

The standard parametrization of the hadron tensor W^ u of deep inelastic scattering (neglect- 
ing the contribution of weak currents) has the form [[/]: 



where p, s, m are the four momenta of the incoming parton (quark), its polarization vector 
and its mass, resp., and q the lepton momentum transfer (— q 2 = Q 2 ). The structure func- 
tions Fi, F 2 , gi, and g 2 depend upon Q 2 and x = Q 2 /2pq. 





+^e^ ua(3 q a s f) —gi + le^pq 01 ^ - p P —) — 
pq pq pq 



(2.1) 
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The structure functions in (1) are energy discontinuities of scattering amplitudes for the 
elastic Compton scattering of a virtual photon off the proton. In analogy with (1) we write: 



with 



= i J d 4 xe tqx < N\T(J^x)J„(0)) \N > 

= {-g^ + -^)Ti + fa-qp-z){pv-qu-~)T2 
qz gz gz 



+te^ a pq a s p —T 3 + it^ aP q a {s p {pq) - /( S g))-T 4 (2.2) 
pq pq 



F x = ~ImT x F 2 = -^-ImT 2 (2.3) 

Z7T 111 



9i = -^ImT 3 g 2 = -^ImT^ (2.4) 

It is easy to see that these amplitudes satisfy the following crossing symmetry relations: 
T\,T 2 and T4 are symmetric with respect to the replacement s —>■ — s (i.e. s — > u channel), 
while the amplitude T3 changes sign. Indeed, the tensor T^ v is symmetric under the inter- 
change of fi and v and q — > —q. But q — ► —q means x — > —x, i.e. s — * u. Since the tensors 
in front of T% and T 2 are symmetric, T\ and T 2 must be even under s — >• u. Correspond- 
ingly, T 3 must be odd, and T 4 whose tensor has an additional power of q, is even again. In 
other words, Ti,T 2 ,T 4 are the amplitudes of the positive signature, and T 3 has the negative 
signature. [] In the following we shall concentrate on this odd signature case T 3 which gives g x . 

One can check the main implications of the signature assignment by simply calculating 
the amplitudes in lowest order perturbation theory. In the Born approximation they take 
the form: 

(2.5) 

TP = e 2 Q ^f- \_ Q2 + ie +_ s _q 2 + J > ( 2 - 6 ) 




T 3 0) - e '( s _Q2 + , e _ s _g 2 + J> ( 2 - 7 ) 



T 4 (0) = (2.8) 

(e q is the electric charge of the initial quark of flavour q.) We have deliberately chosen this 
particular way of writing eqs.(5-8) in order to stress that T\ and T 2 are symmetric under the 

4 All the other kinematical factors (like pq) which sometimes enter into the definition of the structure 
functions do not change the signature (i.e. crossing symmetry) properties of a given amplitude. 
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replacement s 



—s, and T3 is antisymmetric. 



Now let us consider the radiative corrections to the Tj. In the one-loop approximation 
there are two Feynman graphs yielding DL-contributions to (see Fig. 2). However, since 
ultimately we are interested in the discontinuity of the amplitude, we have to consider, be- 
sides the DL-terms, also the «7r-parts of the energy logarithms. When s > 0, only the graph 
2a yields the ^-contribution. The result is: 



n 



(i) 



1 16tt 2 



In' 



■■• /, r '_-2//r^ 
fi z /i z [i A 



S s 2 S 



(i) 



' q 167T 2 



= 2xT{ 
In 2 



(i) 



T 



(1) 



-In 2 



(2.9) 

(2.10) 
(2.11) 



- , - J (2.12) 

where Cp = (N 2 — 1)/2N for the colour group SU(N), g is the QCD-coupling constant 
which we will keep fixed within the DLA, and fi is an infrared mass scale for the transverse 
momentum: 



u< k t 



(2.13) 



Eqs.(9-ll) show that the one-loop expressions, T { , for the functions Tj still have the same 
signatures as the Born terms. As a result, in T\ and T2 the double logs add up, whereas in 
T3 only a term ~ mln(s / fi 2 ) survives. For the energy discontinuity, in both T\ and T3 the 
in pieces in ln 2 (s/ fi 2 ) give the relevant contribution. This pattern will not change if we in- 
clude more rungs in higher orders in g 2 , i.e. when building the usual sum of ladder diagrams. 



However, non-ladder Feynman diagrams can also contribute to Tj in the DLA, and their 
role is quite different for the even- and odd-signature amplitudes. For example, from the 
case of the elastic scattering of quarks it is known that non-ladder graphs contribute 
differently to positive and negative signatured amplitudes. In our case, the non-ladder dia- 
grams first appear (within the DLA) in the two loop approximation, i.e. in the order g 4 of 
perturbation theory. They are illustrated in Fig.3. Let us study these contributions in more 
detail. In particular we want to demonstrate and explain that they do contribute to T 3 but 
not to Ti or T2. 

Following the line of arguments introduced in 0, each of the non-ladder diagrams can 
be viewed as being obtained from a ladder graph by adding non-ladder virtual gluons to it. 
In order to give the double logarithm, these gluons couple only to the side lines of the ladder 
(i.e. not to the rungs), and they must be softer (i.e. have smaller k t ) than all those ladder 
gluon rungs which are embraced by the non-ladder gluon. Otherwise they would destroy 
double logs of the ladder. Morover, those propagators of the ladder to which the non-ladder 
gluon is attached, must be close to the mass shell - its virtuality should be smaller than the 
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transvers momentum square of the non-ladder gluon. In the following, these non-ladder type 
gluons will be called bremsstrahlung gluons. The important consequence of their "softness" 
is that, in the numerator of the expession for the non-ladder graph, the momenta of the 
bremsstrahlung gluons can be neglected compared to the ladder gluon momenta. Further- 
more, the momenta along the side lines of the ladder, are small in comparison with the initial 
momenta q and p. In Fig. 3 the momentum of the bremsstrahlung gluon is denoted by k<i- 
The remaining momenta in all these numerators give, for each diagram, the factor ±4k 2 -pk\ 
(where the sign plus stands for the diagrams in Fig.3a,b, while the minus sign belongs to 
Figs. 3c, d). In front of this numerator we have a tensor which is built from p^, or g^ v . 
These tensors are identical for the graphs (a,b,c,d), and those of the (a',b',c',d')-diagrams 
differ by the replacement \i v and q — > —q. 

The way in which these double logs sum up or cancel, depends upon color and signature. 
First it is easy to see that for a color singlet t-channel different bremstrahlung gluons can- 
cel: for example, in Fig. 3a and c the bremsstrahlung gluons are emitted from the lower left 
external line, and we are summing over different end points. Since we have color zero in the 
t-channel, both contributions cancel. As a result, in this order all double log bremsstrahlung 
contributions in Fig. 3 cancel (this will change when we go to higher order, as we shall discuss 
further below). 

Signature becomes crucial when we consider the vk terms which are essential for the s- 
discontinuity. To be more precise, starting from the order g 4 where the nonladder graphs 
appear, even and odd signature amplitudes start to behave quite differently. For example, 
the graphs in Figs. 3a' and b' have no s-cuts, but those in c' and d' can be redrawn as shown 
in Fig.4, and therefore they contribute to the s- discontinuity. Combining Figs.3a,b with c',d' 
one finds cancellation for even signature, non-cancellation for odd signature. As a result, 
the non-ladder graphs contribute to T3, but not T\ or T%. This is in accordance with the 
conclusion of || for the quark scattering amplitude. 

When we move on to higher orders in g 2 , most of these features can easily be generalized. For 
the remainder of this section we remove the photon lines at the upper end, and we restrict 
ourselves to quark scattering. First let us return to the bremstrahlungs gluons. Starting 
from a fermion ladder with a few gluon rungs, each non-ladder gluon will comprise a subset 
of the rungs. In order to give the maximal number of double logarithms, this bremsstahlung 
gluon must be softer (i.e. have smaller k t ) than any line of this subset of rungs; at the same 
time, the fermion lines to which the bremsstrahlung gluon couples must be softer than the 
gluon itself. This implies that the bremsstrahlung gluons form a "nested" structure. The 
summation of all these graphs is done as described in ||. Consider a typical diagram; pick 
the internal line with the lowest k t . If it is a fermion line, the amplitude can be decomposed 
as indicated in Fig. 5a. If, on the other hand, it is a gluon line, it must be a bremstrahlungs 
gluon which - according to Gribov's theorem - |S| must go from one external leg to another 
(Fig. 5b). The sum of all diagrams, therefore, can be written as shown in Fig. 5c; the first 
term on the rhs denotes the Born term. If we take the derivative with respect to the infrared 
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cutoff /i 2 (which appears only in the line with the smallest momentum), then we arive at the 
infrared evolution equation illustrated in Fig. 5c (whether or not the Born term contributes 
to this derivative has to be decided from case to case). This is the infrared evolution equation 
which has proven to be an efficient instrument for the investigation of the elastic scattering 
of quarks ||, and also in inelastic reactions at high energies |T0| . 

The presence of the second term on the rhs has an important consequence. Even we re- 
strict ourselves to amplitudes with color zero in the t-channel, the colored gluon forces the 
four quark amplitude to be in a color octet. Therefore, this second part appears as an 
inhomogeneous term in the nonlinear evolution equation, and we have to solve a separate 
evolution equation for the color octet amplitude. Furthermore, in the case of a colour sin- 
glet amplitude this term contributes only to the odd signature, as it can be seen as follows. 
Consider, for example, the two diagrams shown in Fig. 6a and b. For the odd-signature 
amplitude, they both come with the opposite sign. Redraw the first diagram as illustrated 
in Fig. 6c, and note that the upper line to which the soft gluon couples has the opposite 
direction: this gives an additional minus sign from the color factor, and the two diagrams 
sum up. Conversely, in the even signature case they cancel. As a result, the nonplanar 
diagrams appear only in the odd signature amplitude. 

Finally, the infrared evolution equation for the (even signature) octet amplitude. Following 
the same arguments as before, we are led to the nonlinear evolution equation whose structure 
we illustrate in Fig. 7. In the second term, only the even signature color octet contributes: 
the gluon has negative signature, and for the subamplitude the even signature configuration 
gives the maximal number of double logs. Combining this with the color factors, one finds 
that only the color octet survives. 

3 Infrared Evolution Equation for g\ 

In this section we apply the methods outlined in the previous section in order to construct 
and solve the infrared evolution equation for the nonsinglet contribution to g%. We will take 
into account both the DL-contributions and the m contributions. 

Let us consider, again, the forward elastic photon-quark scattering amplitude T^ u in DLA. 
As before, q and p are the (collinear) four momenta of the photon and the quark, resp., and 
we are in the regime where s ^> — q 2 = Q 2 ^> p 2 . The structure of the infrared evolution 
equation for quark scattering has been discussed in the previous section, and we now apply 
the same ideas to photon quark scattering. The equation is illustrated in Fig. 8, and its 
origin is easily understood. The upper blob on the rhs denotes the photon-quark scattering 
amplitude with the lower (incoming) legs having transverse momentum square equal to the 
infrared cutoff fi 2 . The lower blob on the rhs represents a quark quark elastic scattering 
amplitude with k\ = fi 2 for all external momenta. In order to translate this into an explicit 
equation, we note that T^ v depends upon the two variables z = ln(s/^ 2 ) and y = ln(Q 2 / fj, 2 ). 
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Therefore, the lhs of Fig. 5 has the two terms: 

BT 8T 8T 

_ n 2 V v — ^ V _|_ ^ V ]\ 

dfi 2 dz dy 

The convolution on the rhs becomes simpler if we write our amplitudes in the Mellin repre- 
sentation w.r.t. s (z = ln(s//x 2 )): 

/ico Ai,\ q 
fi( r^ Rl (uj,y), (3.2) 

where u denotes the angular momentum j, and £ is the signature factor: 



e l " w + r, 



ew = 2 " (3.3) 

with r 3 = — 1. As usual, the integration contour runs to the right of all singularities in the 
a;-plane. A similar representation is used for the quark-quark scattering amplitude, and the 
transform is denoted by /q . The infrared evolution equation of Fig. 8 then takes the form: 

+ = ^^/o (+) (3-4) 
ay 8tt 2 

for i—1,2 and 

f)R 1 

for the amplitude T 3 . Eq.(3.4) has been discussed and solved in 0. 

In order to solve eq.(3.5) for R^, we need to know /g . This is the quark scattering amplitude 
discussed in the previous section. It satisfies the evolution equation of Fig. 5 [§]: 

( _) N 2 -l g 2 N 2 - 1 g 2 1 (+) 1 /,(-)/ ,\ 2 ,„ R . 

The second term on the rhs is due to the signature changing contributions of gluon bremsstrahlung 
which lead us to define an even-signature quark quark scattering amplitude /g (u>) with color 
octet quantum number in the t-channel. This amplitude has the infrared evolution equation 
shown in Fig. 7: 

Its solution follows from the discussion given in ||: using the transformation 

ft\u)^Ng 2 ^-\nu{u) (3.8) 
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the Riccati equation (3.7) is equivalent to the linear differential equation 

du 2 du 1 



, , r-rM = (3.9) 

dz 2 dz 2N 2 K ' 



where 



z = u/u , ujq = ^Ng 2 /%n 2 . (3.10) 



After a trivial transformation this equation is solved by a parabolic cylinder function. As a 
result, /g has the form: 

/ 8 (+) H = Ng 2 ^ln(e^D p (z)) 

with 

p=-l/2N 2 . (3.12) 
With this we return to (3.6) and obtain for /q ^: 



(the minus sign in front of the square root follows from the requirement that, for large u, 
the solution has to match the Born approximation). 

Having found the function /q \uj), we are now able to solve eq.(3.5): 

R 3 (aj, y) = C(uj)e^ + ^ ] ^l^ y (3.14) 

with some unknown function C(oS). This function is determined if we impose the matching 
condition [0] : 

T 3 (z,y = 0)=f 3 (z) (3.15) 

where T 3 (z) is the corresponding invariant amplitude with a (nearly) on-shell photon: — q 2 = 
fi 2 . Since this amplitude no longer depends upon Q 2 , its infrared evolution equation (which 
is the analogue of eq.(3.5))) has no y- derivative. On the other hand, its Born term now has 
a dependence upon /i 2 which is no longer negligeable and leads to a nonzero contribution on 
the rhs (Fig.8). Denoting the Mellin transform of T 3 by R 3 the matching condition (3.15) 
simply becomes C(uj) = R 3 , and the equation for R 3 reads: 

ooR 3 = c 3 + -^R 3 ft ) . (3.16) 



The boundary condition for T3 is: 

Tz(z = 0) = T 3 Born = C3 

R 3 = - (3.17) 

UJ 

where T 3Born is given in (2.7) and leads to c 3 = 2e 2 . 
Inserting the solution for R 3 = C(u) into (3.15) we obtain: 

^--^rs^^w"^ <3 - i8) 

In order to obtain our final result, g±, we have to take the discontinuity in s. With the 
signature factor = vkuj/2 and the variables x = Q 2 /s and Q 2 we arrive at 

q 2 ) = fl r « (3 , 9) 

In the limit where 1/x is much larger than Q 2 //! 2 , the leading contribution comes from the 
rightmost singularity in the tu-plane. 

4 Discussion 

The most interesting aspect of our result (3.19) is the small-x behaviour. The leading 
singularity in the u;-plane is the branch point due to the vanishing of the square root in /q _ ' ) 
in (3.13). As a first approximation, let us take the number of colours iV to be large. As it 
has been discussed in |§, /g can then be approximated by the first term on the rhs of 
(3.7) (i.e. its Born term), and the location of the branch point follows from the condition: 

q 2 N q 2 

= 1--^—- 1 + — 1 4.1) 

(note that our u = j, i.e. it differs from the u = j — 1 which is commonly used in the context 
of the BFKL Pomeron and the unpolarized gluon structure function). Expanding in inverse 
powers of N we find for the first two terms: 



UJ — u/ ^ = UJ 



1 



2N 2 

where 



+ (4-2) 



„W = (4.3) 

V 4iV7T 2 ^ ^ 
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is the rightmost singularity for the flavor nonsinglet part of Fx found in 0. Lut us stress 
that the leading singularity of the negative signature amplitude (gx) is shifted to the right 
compared to the the positive signature one M. Fistly such a behaviour of the negative 



signature amplitudes were noticed in [11] for the elastic scattering in QED 



A more accurate determination of the branch point singularity leads, instead of (4.2), to: 

W H ~ CJ (+) . 1.04 

w 0.41 (4.4) 

for a s = 0.18. As a result, the power of 1/x of the nonsinglet piece of g± is stronger than 
that of F\ by about 4 %. This enhancement has its origin in the non-ladder graphs which 
where absent in the calculation of Fx- In the HERA region, the enhancement of gx relative 
to Fx [0 is approximately 

gx/Fx « 1. W +) 0.04[ln(l/x) + ^ln(Q7/i 2 )] 

« 1.13 (4.5) 

for x = 10 -3 , Q 2 = 20GeV 2 , and /i 2 = AGeV 2 . This estimate may serve as a crude method 
of estimating the non-singlet g™' s ' in the small-x region, once FJ 1 ' 8 ' has been measured. 

Let us finally see how our result (3.19) is related to the GLAP expression. To this end 
we take Q 2 /fi 2 to be much larger than 1/x; the exponent /q \uj)/8tt 2 plays the role of the 

2 • 



anomalous dimension ^y qq . For y/aT s <C uj <C 1 (i.e. g /to is small) we obtain from (3.13): 

. . oc s N 2 — 1 , t , 

^ = 2 *£-2fT (4 - 6) 

which agrees with the singular part of the quark anomalous dimension. In a straightforward 
extrapolation of the GLAP approximation into the small-x region ]T2"| , |T3| , 13 one would use 



this singular term in the anomalous dimension and performe a saddle point analysis. This 



has lead to the conclusion that at small x gx ~ exp(y const ln(Q 2 /A 2 ) ln(l/x)), in analogy 
with the small-x behavior of the gluon structure function (at fixed a s ). In contrast to this, 
our result in (3.19) gives a quite different answer for the uj dependence of the quark anoma- 
lous dimension and hence for the small-x behavior of gx- For a numerical estimate of the 
difference we simply combine (4.5) with the result for Fx which was obtained in || (note 
that for the flavor nonsinglet case the 7 gg is the same for the polarized and the unpolarized 
structure function): for the non-singlet quark structure function with the initial condition 
q(x, Qq) oc 5(1 — x) f\ a comparison has been made of the GLAP-evolution at small-x and 
the double logarithmic formula derived in 0, and a difference of up to a factor of 10 (at 
typical HERA values for x and Q 2 ) was found. Together with (4.5), this leads to a slightly 
stronger discrepancy in g\\ a factor of ten or even more. 



5 It does not seem very plausible to assume that gi(x,Q^) is singular at x — > 0, therefore we think that 
the somewhat oversimplified (5-function ansatz is justified. 
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An even stronger discrepancy exists between our result and the Regge prediction [16], [T7 |. 
For the flavor nonsinglet (isotriplet) exchange the small-x behavior is given by a\ exchange. 
It has the form (l/x) aR , and the a\ intercept an is believed to lie somewhere between —0.5 
and [111], i.e. the flavor nonsinglet part of g\ is predicted to be regular as x — > 0! 



Having found that the double logarithms in the small-x region are numerically important, 
we have to adress the slightly more involved analysis of the flavor singlet case. Here we are 
facing the mixing between gluons and quarks: as it can be seen from 0] or, more directly, 



from p2|, |17| , the small-x behavior of the polarized gluon structure function is less singu- 
lar than the unpolarized one and, therefore, competes with the quarks. This analysis is in 
progress, and results will be presented in a forthcoming paper. 
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Figure Captions: 



Fig.l: The Born approximation for T^ v . 
Fig. 2: The one-loop approximation for T^ v . 
Fig. 3: Lowest order non-ladder gluons in T^. 

Fig. 4: The right hand cut energy discontinuities of the diagrams Figs. 3a, b, and c': in 
the even signature amplitude, a and c' cancel (similary a and d', which is not shown in 
Fig. 4), whereas in the odd signature case they are coming with the same sign and add up. 

Fig. 5: The line with lowest k t in a general non-ladder diagram for the quark quark scat- 
tering amplitude, (a) a quark line, (b) a gluon line, (c) the infrared evolution equation for 
the sum of all non-ladder diagrams. The sum in front of the second term on the rhs denotes 
the different possibilities of attaching bremsstrahlung gluons to the external lines. 

Fig. 6: Cancellation of bremstrahlung gluons in the even signature amplitude. 

Fig. 7: The infrared evolution equation for the octet even signature quark quark scattering 
amplitude. The sum in front of the second term on the rhs denotes the different possibilites 
ways of attaching bremsstrahlungs gluons to the external lines. 

Fig. 8: The infrared evolution equation for the photon quark scattering amplitude. The 
lower blob on the rhs denotes the quark quark scattering amplitude from Fig. 5c. 
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